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Abstract 

Recent results on the energy spectra in QCD jets are reported. Within the Modified Lead- 
ing Logarithmic Approximation (MLLA) and the Local Parton Hadron Duality (LPHD) 
model one finds a very good description of the e + e~ data from the lowest up to the LEP 
1.5 energies. A model with fixed a s can be excluded, already in certain finite energy 
intervals. The fits also extrapolate smoothly into the region of small particle energies, in 
particular, the data follow a scaling prediction for the low energy limit derived from the 
colour coherence of the soft gluon emission. 
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1 Introduction 



2 Moment Analysis of Particle Spectra 



One of the important predictions of perturbative 
QCD on the intrinsic structure of jets concerns 
the energy spectra of particles. The particles of 
low energy E are not multiplied with increasing 
jet energy Ej et because of colour coherence in the 
cascading process and this yields the bell shaped 
spectrum in the variableJ^p log , the so-called 
"hump-backed plateau" M 

Predictions on the spectrum have been car- 
ried out in the MLLA which takes into account the 
leading double logarithmic (DLA) results and all 
next-to-leading corrections of order */a7. Terms 
of higher order are included as well, although not 
completely, which allow to fulfill the initial con- 
dition for the parton cascade at threshold. These 
predictions, at the parton level, involve only two 
parameters: the QCD scale A which determines 
the running of the coupling a s and the transverse 
momentum cut-off Qq of the gluon emission. Re- 
markably, the observed hadron spectra are rather 
well described by the spectrum of partons if a low 
value Q ~ 250 MeV of the order of hadronic 
masses is used ancLthis observation has led to the 
LPHD hypothesis U. As a justification of this ap- 
proach is not yet available at a fundamental level 
it seems important to determine its range of ap- 
plicability and its limitations. Also one would like 
to know to what extent the predictions are sen- 
sitive to QCD as the underlying theory. In this 
contribution the sensitivity of the energy spectra 
to the running a s and the colour coherence - two 
properties specific to QCD as a field theory 
addressed. More details can be found in rcfs.Li 



To analyse the effect of the running a s it is 
convenient to work with the moments of the £- 
spectrum of particles which evolve independently 
of each other with energy. They are defined by 
< £ 9 (T,A) >= / d££ q D(l,Y,\)/J\f where we use 
the logarithmic variables Y — log(£^ e t/Qo) and 
A = log(<5o/A); the normalization is by the mul- 
tiplicity N. We consider finally the cumulant mo- 
ments K q which are obtained from the < £ 9 > mo- 
ments by Ki =< £ >= £, K 2 = a 2 =< (f - £) 2 >, 
K 3 =< (£ - C) 3 >, K 4 =< (C - 4 > -3a 4 , 
. . . ; also one introduces the reduced cumulants 
n q = Kq/a q , in particular the skewness s = K3 
and the kurtosis k = K4. The cumulant moments 
K q (Y, A) behave at high energies like 



K q (Y,X) = Kq(Y ,X) + 



dy 



du> 



7u< A)] 



(1) 



where Yq is the initial energy. Here "f UJ [a s (y, A)] 
denotes the anomalous dimension of the Laplace 
transform of the ^-spectrum dn/d^ = £>(£,Y,A). 
Eq. (|l|) shows directly the sensitivity of the cumu- 
lant moments to the running of a s . In particular, 
for fixed coupling K q ~ Y at high energies. In 
our applications, we take into account the initial 
condition at threshold for the jet to consist of only 
one parton which implies K q = at Yq — 0. 

The MLLA prediction for the £ 9 moments can 
be written for arbitrary Qq and A asQ: 



k=0 



(2) 



where N r , N 2 , L ( k q) and R k q) arc known functions 
of a = lli\r c /3+2n//3iV c 2 , b = (llN c -2n f )/3, and 
A where n f , Nc denote the numbers of flavours 
and colours. Comparing these-formulae to the ex- 
perimental data on moments, u best agreement is 
found for the "limiting spectrum" , where the two 
parameters coincide, i.e., Qo — A, or A = 0. In this 
case the formulae simplify and the moments can 
be expressed! in terms of the parameter B = a/b 
and the variable z = y/l6N c Y/b as: 

' ° ' = P Q iq) (B + l } B + 2,z) + 



2lB+2{z) -pW(B + l,B + 2,z) 



z I B+ \{z)' 



(3) 



where and P[ q ^ are polynomials of order 
2(q~ 1) in z. 

These expressions extrapolate smoothly to 
threshold, where they are determined by the ini- 
tial condition for a single parton. Similarly one 
can. derive the complete MLLA results for fixed 
a s u, which simplify for high energies to 
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where r] = a"f^/8N c , 7o = V7o + V 2 and 7o = 
y/2N c a a /n. 

3 Relating parton and hadron spectra 

Strictly speaking the QCD results are only reliable 
for E Qq. In order to determine the moments 
one has to integrate the spectra over the full range 
and to extrapolate to the soft limit. The experi- 
mental particle densities are usually given as func- 
tion of momenta. The same kinematic limits of 
parton and hadron spectra is obtained if the same 
effective mass Qo is assigned to the charged par- 
ticles which also limits the energy of the partons 
(Eh = E p > Q ). 

Furthermore, a common behaviour of the par- 
ton and hadron spectra near the boundary can be 
obtained if the spectra are related bylM 

^ = ^ sW)A) ( 6 ) 



dph 



dp P 



at the same energy E or £e = £ = \og(Ej et /E), 
where E h = y/p\ + Ql and E p = p p ; K h de- 
notes a normalization parameter. With this choice 
both spectra vanish linearly in £; this corresponds 
to the invariant distribution of hadrons E-jjsrj: ap- 
proaching a constant value in the soft limit p —> 0, 
as suggested by the experimental data. The pre- 
scription ^ represents a minimal extension of the 
QCD prediction for E ^> Qo towards the bound- 
ary at E = Qo but is not unique; other prescrip- 
tions using a smaller value of Qo ~ m-K have also 
been applied (see discussion inQ). 

4 Comparison with experimental data 

The moments < £ 9 > are determined i from the 
spectra Edn/dp vs. after appropriate trans- 
formation of the measured x p = 2p/^/s spectra 
of charged particles using E = ^/p 2 + Qq and 
therefore depend on Qo- First a fit of the two 
parameters Qo and A is obtained by comparing 
the moments for a selected Qo with the theoreti- 
cal predictions from Eq.(^) for different A. The 
best agreement with the data was obtained for 
Qo fa A » (270 ± 20) MeV. It was not possible 
to obtain a satisfying description for the fixed a s 
case. Choosing 7 = 0.64 and the same Qo = 
270 MeV a good description of the multiplicity 
Af and the slope of £ vs. Y could be obtained but 
not for the other quantities. In Fig. 1 we show the 
evolution of the first cumulant moments K q with 
enery for running a s and fixed a s . Note that these 
predictions depend only on the two parameters Qo 
and A which actually coincide in the fit. The ab- 
solute normalization of the moments is given at 
threshold Y o =0 by K q =0. 

As this is an application of perturbative QCD 
to an extreme limit one could make the weaker 
(more conventional) assumption and choose a 
higher starting energy, say Yo — 2 (y^s ss 4 GeV). 
In this case for each moment an extra free param- 
eter K q (Yo) had to be introduced which, however, 
would not improve the fit for running a s essen- 
tially; a backward evolution from Yo would again 
reproduce approximately the initial condition at 
threshold. Therefore, the initial condition yields 
the highly constrained fit with only two parame- 
ters (the normalization parameter only enters 
the multiplicity AA, not shown here) . 

In case of fixed a s one may ask the same ques- 




1 - 



0.2 - 

^ 





012345678 
Y 



Figure 1: The first four cumulant moments K q of the 
charged particles' energy spectra Edn/dp vs. fg, shown 
as a function of Y for Qo = 270 MeV. Predictions of the 
Limiting Spectrum of MLLA with running a B (solid lines), 
of MLLA with fixed a s (= 0.214) (dashed lines) and of 
MLLA with fixed a s normalized by hand to the experimen- 
tal point at \fs = 44 GeV (dotted lines) are also shown; in 
all cases rij = 3 (from Ref . 4) . 
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Figure 2: Charged particle inclusive momentum distri- 
bution at LEP-1.5 from ALEPH (diamonds), DELPHI 
(squares) and OPAL Collaborations (triangles) in compar- 
ison with theoretical predictions of the Limiting Spectrum 
with Qo = 270 MeV (solid line). Dashed lines show the 
predictions of the Limiting Spectrum after correction for 
kinematical effects (from Ref. 5). 

tion, whether the fit could be improved by choos- 
ing a higher starting value Yq, i.e. to assume a 
transition from a running coupling to a fixed cou- 
pling regime at Yq (This would correspond to a 
parallel shift of the dashed curve to the data at 
Yq). As can be seen in Fig. 1 the energy depen- 
dence of the fixed and running a s curves is quite 
different so that an improvement cannot be ob- 
tained in a finite energy interval. Rather, such 
a procedure can be excluded for a Y interval of 
about 1-2 units. An example is shown by the dot- 
ted curve which corresponds to Yq = 4.3 (y^s = 
44 GeV) . The trend at yet larger energies is again 
different, in particular for the higher moments. 

The dependence on the. number of flavours rif 
has been studied as wellurcl. As the jet evolution 
is dominated at all energies by the gluon emission 
at the lower transverse momenta, the choice nt = 
3 is a good approximation up to LEP energies. 

The recent LEP-1.5 data are also well ac- 
counted for la by the same theoretical scheme as 
is shown in Fig. 2. Here the experimental £- 
spectrum is shown in comparison with the lim- 
iting spectrum (setting £ p = which terminates 
at Y ~ 5.3 and after rescaling £e — > £ p with re- 
lation (^|) which takes into account the boundary 
effects. Also the moments of the distribution and 
the position of the maximum are well reproduced 
by the theoretical predictions at this energy. 
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Figure 3: Invariant density Edn/d 3 p of charged particles 
in e+e - annihilation as a function of the particle energy 
E = y/ P 2 + Qq at Qo = 270 MeV at various cms energies 
in comparison with the MLLA predictions (from Ref. 6). 



5 Soft limit of energy spectrum 

Having observed the good fit of the QCD-LPHD 
prediction not only towards high but also towards 
low cms energies yfs we now turn to the behaviour 
of the fits at low particle energy E. As the emis- 
sion of very soft gluons from all other partons in 
the jet is coherent the production rate of such par- 
ticles is nearly energy independent a. Indeed, the 
analytic results in the DLA and MLLA show the 
scaling behaviour of the spectrum in the soft limit 
and the production rate is essentially determined 
by the colour charge of the initial partons U. Re- 
markably, the observed hadrons follow the trend 
of this prediction as can be seen in Fig. 3 which 
shows the invariant particle density as a func- 
tion of particle energy E (assuming the effective 
charged particle mass Qq = 270 MeV). The data 
from all energies in e + e~ annihilation approach a 
common limit for E — > Qq (p — > 0). The curves 
which represent the MLLA calculation show again 
the scaling limit and reproduce the dependence 
on energy E and cms-energy yfs rather well. This 
suggests that LPHD persists towards very low par- 
ticle energies. Various further tests of this picture 
have been proposedQ by studying the dependence 
of the soft limit Jo on the colour charge of the 
primary partons in e + e~ — > 3 jets, deep inelastic 
scattering processes and hadronic collisions. 



6 Conclusions 

The particle energy spectra are well described by 
the QCD parton shower in the MLLA assuming a 
close similarity between partons and hadrons ac- 
cording to the LPHD approach. This similarity 
appears to work not only for the highest but also 
down to rather low cms energies as well as to low 
particle energies. This agreement is by no means 
trivial. For example, replacing the running a s by 
the fixed a s yields predictions quite incompatible 
with experiment, also if the threshold for the on- 
set of a fixed a s regime is increased. The scaling 
prediction derived from colour coherence for the 
soft particles works well for hadrons. 

This may suggest □ to view the final stage of 
hadron production through resonances, as explic- 
itly incorporated into Monte Carlo hadronization 
models, as being dual to the partonic evolution 
with running coupling down to the low scale of a 
few hundred MeV. 
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